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Abstract 



According to Feynman, the universe consists of two parts - the system 
in which we are interested and the rest of the universe which our measure- 
ment process does not reach. Feynman then formulates the density matrix in 
terms of the observable world and the rest of the universe. It is shown that 
coupled harmonic oscillators can serve as an illustrative example for Feyn- 
man's "rest of the universe." It is pointed out that this simple example has 
far-reaching consequences in many branches of physics, including statistical 
mechanics, measurement theory, information theory, thermo-field dynamics, 
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quantum optics, and relativistic quantum mechanics. It is shown that our 
ignorance of the rest of the universe increases the uncertainty and entropy in 
the system in which we are interested. 

I. INTRODUCTION 

Because of its mathematical simplicity, the harmonic oscillator provides soluble models 
in many branches of physics. It often gives a clear illustration of abstract ideas. In his book 
on statistical mechanics |l| , Feynman makes the following statement on the density matrix. 
When we solve a quantum-mechanical problem, what we really do is divide the universe into 
two parts - the system in which we are interested and the rest of the universe. We then 
usually act as if the system in which we are interested comprised the entire universe. To 
motivate the use of density matrices, let us see what happens when we include the part of 
the universe outside the system. 

The purpose of this paper is to study Feynman's rest of the universe and related problems 
of current interest using a pair of coupled oscillators. Starting from the classical mechanics 
of harmonic oscillators, we formulate the symmetry of the oscillator system in terms of the 
group Sp(2) which is of current interest and which is locally isomorphic to the (2 + 1)- 
dimensional Lorentz group. This symmetry is then extended to the quantum mechanics of 
coupled oscillators. This allows us to study measurable and unmeasurable variables in terms 
of the two oscillator coordinates. The unmeasurable variable constitutes Feynman's rest of 
the universe. We shall study the effects of this unmeasurable variable to the measurable 
variable in terms of the uncertainty relation and entropy. 

In Sec. we reformulate the classical mechanics of two coupled oscillators in terms 
of the Sp(2) group. The symmetry operations include rotations and squeezes in the two- 
dimensional coordinate system of two oscillator coordinates. In Sec. fTJ, this symmetry 
property is extended to the quantum mechanics of the coupled oscillators. In Sec. [TV], 
we use the Wigner phase-space distribution function to see the effect of the unobservable 
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variable on the uncertainty relation in the observable world. In Sec. [V], we use the density 
matrix to study the entropy of the system due to our ignorance of the unobservable variable. 
In Sec. |VT], it is shown that the system of two coupled oscillators can serve as an analog 
computer for many of the physical theories and models of current interest. Section |V11| 
contains some concluding remarks. 



II. COUPLED OSCILLATORS IN CLASSICAL MECHANICS 

Two coupled harmonic oscillators serve many different purposes in physics. It is widely 
believed that this oscillator problem can be formulated into a problem of quadratic equation 
of two variables, and the quadratic equation can be separated by a simple rotation. It is 
true that the problem can be reduced to a quadratic equation, but it is not true that this 
equation can be solved by one rotation. Indeed, in order to understand fully this simple 
problem, we have to employ the SL(4,r) group which is isomorphic to the Lorentz group 
0(3,3) with fifteen parameters 0. 

In this paper, we do not need all the symmetries available from the 0(3,3) group. We 
shall use one of the 0(2,l)-like subgroups which is by now a standard language in physics 
in all branches of physics. Let us consider a system of coupled oscillators. The Hamiltonian 
for this system is 

H = - \— p\ + — pi + Ax\ + Bx\ + Cx x x 2 \ . (2.1) 
2 Imi 1712 J 

By making scale changes of X\ and x 2 to (m 2 /m 1 ) 1 / 4 Xi and (m 1 /m 2 ) 1 / 4: x 2 respectively, it is 

possible to write the above Hamiltonian in the form [[|[|] 

H = hi + \ { Ax * + Bx " + CxiX2 } ' (2,2) 

with m = (m 1 m 2 ) 1 ^ 2 . We can decouple this Hamiltonian by making the coordinate trans- 
formation: 

yi \ /cos (a/2) — sin(a/2)\ / x\ \ 
y 2 J \sin(a/2) cos(a/2) / \x 2 ) 
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Under this rotation, the kinetic energy portion of the Hamiltonian in Eq. (|2.3|) remains 
invariant. Thus we can achieve the decoupling by diagonalizing the potential energy. Indeed, 
the system becomes diagonal if the angle a becomes 

C 

tana = -^— -. (2.4) 

This diagonalization procedure is well known. What is new in this note is to introduce the 
new parameters K and rj defined as 



K = JAB - C 2 /4, 



A + B + J(A-B) 2 + C 2 



In terms of this new set of variables, the Hamiltonian can be written as 

h =^{p*+pI} + Y + e "^} ' (2 - 6) 

with 



a a 
yi = xicos- -x 2 sm-, 

y 2 = xism- +x 2 cos-. (2.7) 

In this way, we can study the symmetry properties of the coupled oscillators systematically 
using group theoretical methods. The coordinate rotation of Eq. (Ol) is generated by 
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The rj variable changes the scale of y\ in one direction while changing that of y 2 in the 
opposite direction. If y 1 is expanded then y 2 becomes contracted so as to preserve the 
product ?/i?/2 • This is called the squeeze transformation. The squeeze operation which 
changes the scales of y\ and y 2 is generated by 

* = -5 (2 ' 9) 

From the linear transformation of Eq. (|2.3| ), this expression can be written as 



= S\ cos a — S 2 sin a, 



(2.10) 



with 

Si = -- L— - X2 — \ 

2 I dxi dx 2 J ' 

Indeed, the generators J , Si and 52 satisfy the commutation relations: 

[Jo, Si] = *S 2 , [Jo, S 2 ] = -iS 1; [Si, S 2 ] = zJ„, (2.12) 

This set of commutation relations is identical to that for the group SU(1,1) which is locally 
isomorphic to the (2 + l)-dimensional Lorentz group Q. If the problem is extended to 
the four- dimensional phase space consisting of the xx, £2, Pi, and p 2 variables, the symmetry 
group is Sp(4) which is locally isomorphic to 0(5,2) [f§-f|. These groups are known to provide 
the standard language for the one- mode and two-mode squeezed states respectively [§]], in 
addition to their traditional roles in other branches of physics, including classical mechanics, 
nuclear, elementary particle and condensed matter physics. 

III. QUANTUM MECHANICS OF COUPLED OSCILLATORS 

If yx and y 2 are measured in units of (mK) 1 ^, the ground-state wave function of this 
oscillator system is 




The wave function is separable in the yx and y 2 variables. However, for the variables x\ and 
x 2 , the story is quite different. 

The key question is how the measurement or non-measurement of one variable affects 
the world of the other variable. If we are not able to make any measurement in the x 2 
space, how does this affect the quantum mechanics in the x\ space. This effect is not trivial. 
Indeed, the x 2 space in this case corresponds to Feynman's rest of the universe, if we only 
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know how to do quantum mechanics in the x\ space. We shall discuss in this paper how we 
can carry out a quantitative analysis of Feynman's rest of the universe. 
Let us write the wave function of Eq.(|3.1|) in terms of x\ and x 2 , then 



ip v {x 1 ,x 2 ) = -^exp j-i 



e '(Xi cos x 2 sin — J 



-e v (xi sin — + x 2 cos ^] 



(3.2) 



If rj = 0, this wave function becomes 



\/7r 



2 ( x i + X 2> 



(3.3) 

For other values of 77, the wave function of Eq. (|3.2| ) can be obtained from the above expression 
by a unitary transformation generated by the operators given in Eq. (|2.8|) and Eq. (|2.11|) . We 
should then be able to write Eq.( [3.2|) as 



^m im2 (a,?7)0 mi (a;i)0 m2 (x2), 



(3.4) 



where 4> m (x) is the m th excited-state oscillator wave function. The coefficients A mim2 (rj) 
satisfy the unitarity condition 



J2 l^m im2 (a,^)| 2 = 1. 



(3.5) 



It is possible to carry out a similar expansion in the case of excited states. 

The question then is what lessons we can learn from the situation in which we are not 
able to make measurements on the x 2 variable. In order to study this problem, we use the 
density matrix and the Wigner phase-space distribution function. 



IV. WIGNER FUNCTION AND UNCERTAINTY RELATION 

In his book, Feynman raises the issue of the rest of the universe in connection with the 
density matrix. Indeed, the density matrix plays the essential role when we are not able to 
measure all the variables in quantum mechanics in the present case, we assume that 
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we are not able to measure the x 2 coordinate. It is often more convenient to use the Wigner 
phase-space distribution function to study the density matrix, especially when we want to 
study the uncertainty products in detail [1],[4|] . 

For two coordinate variables, the Wigner function is defined as Q 

W(x 1 ,x 2 ;pi,p 2 ) = (^) J exp{-2i(p 1 y 1 +P2Z/2)} 

x i)*(x x +yi,x 2 + 2/2)^(^1 ~yi,X2- y 2 )dyxdy 2 . (4.1) 

The Wigner function corresponding to the wave function of Eq. (|3.2| ) is 

W(xi,x 2 ;px,p 2 ) = Q^j exp |-e J? (xi cos ^ - x 2 sin ^) 2 - e~ v (xi sin ^ + x 2 cos ^) 2 

cx cx cx cx 1 

-e~"(piCOS- -p 2 sin-) 2 - e v (p 1 sin - + p 2 cos -) 2 J . (4.2) 

If we do not make observations in the x 2 p 2 coordinates, the Wigner function becomes 



W(xi,p t ) = J W(x 1 ,x 2 ;pi,p 2 )dx 2 dp 2 . (4.3) 
The evaluation of the integral leads to 



, 1/2 

W(xi,x 2 \p 1 ,p 2 ) 



1 



[ 7r 2 (l + sinh 2 7] sin 2 a) j 

x exp { - ( — X A + — V A ) ) . (4.4) 

[ y cosh r] — sin r] cos a cosfi rj + sin r\ cos a J J 

This Wigner function gives an elliptic distribution in the phase space of x\ and p\. This 
distribution gives the uncertainty product of 

(Ax) 2 (Ap) 2 = -(1 + sinh 2 77 sin 2 a). (4.5) 

This expression becomes 1/4 if the oscillator system becomes uncoupled with a = 0. Because 
x\ is coupled with x 2 , our ignorance about the x 2 coordinate, which in this case acts as 
Feynman's rest of the universe, increases the uncertainty in the x\ world which, in Feynman's 
words, is the system in which we are interested. 
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V. DENSITY MATRIX AND ENTROPY 

Since the Wigner function is constructed from the density matrix, it is straightforward 
to show [[§ 

Tr[p(x 1 ,x 1 )] = J W(x 1 ,pi)dx 1 dp 1 , (5.1) 

and 

Tr[p 2 ] = 2n J W 2 (x 1 ,p 1 )dx 1 dp 1 . (5.2) 

If we compute these integrals, Tr(p) = 1, as it should be for all pure and mixed states. On 
the other hand, Tr(p 2 ) becomes 

Tr(p 2 ) = 1/(1 + sinh 2 r? sin 2 a) 1/2 , (5.3) 

which is in general less than one. This gives a measure of impurity and also the degree of 
the effect of our ignorance on the system in which we are interested. 

Let us translate this into the language of the density matrix. If both x± and X2 are 
measured, the density matrix is 

p(x t , x 2 ; x[, x' 2 ) = ip(x u x 2 )ip*{x' 1 , x' 2 ). (5.4) 

In terms of the expansion of the wave function given in Eq. (|3.4j) , 

p(xi,x 2 ;x[,x 2 ) = X Amim 2 (a>,r])A* nin2 (a,r]) 

^ 0mi (%l)4 ) m2 

MKtMKifa)- ( 5 - 5 ) 

If both variables are measured, this is a pure-state density matrix. 

On the other hand, if we do not make a measurement in the x 2 space, we have to construct 
the matrix p(xi,x[) by taking the trace over the x 2 variable: 

PO^H/^,^;,*^, (5.6) 



S 



Then the density matrix p{x\^x' 1 ) takes the form 

p(x 1 ,x[) = CmnO, »7)0 m (a;i)0*(xi), (5.7) 

m,n 

with C mn (a,rj) = A mjfc (a, ??)^4.* fc (a!, 77). The matrix C mn (a,rj) is also called the density 
matrix. The matrix C mn is Hermitian and can therefore be diagonalized. If the diagonal 
elements are p m , the entropy of the system is defined as |7||| 

S = -^p m ln(p m ). (5.8) 

m 

The entropy is zero for a pure state, and increases as the system becomes impure. Like 
Tr(p 2 ), this quantity measures the effect of our ignorance about the rest of the universe. 

It is very important to realize that the above form of entropy can be defined irrespective 
of whether or not the system is in thermal equilibrium. As soon as we have the entropy, we 
are tempted to introduce the temperature. This is not always right. The temperature can 
only be introduced into the system in thermal equilibrium. 



VI. PHYSICAL MODELS 



There are many physical models based on coupled harmonic oscillators, such as the Lee 



model in quantum field theory |9fl, the Bogoliubov transformation in superconductivity [10 



two- mode squeezed states of light ||^,[Tl]], the covariant harmonic oscillator model for the 



parton picture jTJ), and models in molecular physics [13]. There are also models of current 



interest in which one of the variables is not observed, including thermo-field dynamics 14 



two-mode squeezed states [jTJ,!^], the hadronic temperature [17], and the Barnet-Phoenix 



version of information theory [18]. They are indeed the examples of Feynman's rest of 
universe. In all of these cases, the mixing angle a is 90°, and the mathematics becomes 
much simpler. The Wigner function of Eq.(4~2") then becomes 



W(x 1 ,x 2 ;p 1 ,p 2 ) 



71 



exp 



1 r 
2 



e v (xt -x 2 ) 2 + e~ ri {x 1 + x 2 ) 2 

+e- r '(p 1 -p 2 ) 2 + e^p 1 +p 2 ) 2 ]}. 
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(6.1) 



This simple form of the Wigner function serves a starting point for many of the theoretical 
models including some mentioned above. 

If the mixing angle a is 90°, the density matrix also takes a simple form. The wave 
function of Eq.( |3.2| ) becomes 

^(xi,x 2 ) = -^=exp |-i e n (xi - x 2 ) 2 + e" n (x 1 + x 2 ) 2 j . (6.2) 



As was discussed in the literature for several different purposes ||,|(],^l|], this wave function 
can be expanded as 

ip ri {x 1 ,x 2 ) = — \— V(tanh77) fc </>fc(:ri)0fc(:E2)- (6.3) 
cosh 7/ k 

From this wave function, we can construct the pure-state density matrix 

p v (xi,x 2 ; x[, x 2 ) = if)r,{xi, £2)^04, x' 2 ), (6.4) 
which satisfies the condition p 2 = p: 

prj (xi 1 x 2 ] x-y , x 2 ) — I p^ (x\ , x 2 5 x^ , x 2 )prj (^i > % 2 j X\ , x^dxy dx 2 . (6.5) 



If we are not able to make observations on the x 2 , we should take the trace of the p matrix 
with respect to the x 2 variable. Then the resulting density matrix is 

fM*»*,)W*»«)}' (6-6) 

If we complete the integration over the x 2 variable, 

p^x^x'y) = (^ nco ^2r]) ) exp { - i^ Xl +x 'i) 2 / cos H 2r l) + [xi - x[) 2 cosh 2rf\ J . (6.7) 
The diagonal elements of the above density matrix is 

/ l \ 1/2 

P(Xl ' Xl) = \ tt cosh 2 V J 6XP ("^ C ° Sh 2 V ' (6 ' 8) 

With this expression, we can confirm the property of the density matrix: Tr(p) = 1. As for 
the trace of p 2 , we can perform the integration 
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(6.9) 



1 V 



y cosh 77 J ' 

which is less than one for nonzero values of rj. 

The density matrix can also be calculated from the expansion of the wave function given 
in Eq.( |6.3[ ). If we perform the integral of Eq. flS.6p , the result is 

Pv (x 1 ,x[)= (— L- ) ^(tanh^^^O^K), (6.10) 



cosh 77 



k 



which leads to Tr(p) = 1. It is also straight-forward to compute the integral for to Tr(p 2 ) 
The calculation leads to 



Tr(p 2 ) = (l/cosh^) 4 ^(tanhr7) 4fe . (6.11) 



The sum of this series is (1/ cosh rj) 2 , which is the same as the result of Eq.flfT 

This is of course due to the fact that we are not making measurement on the X2 variable. 
The standard way to measure this ignorance is to calculate the entropy defined as [§] 

S = -Tr (p\n(p)) . (6.12) 

If we use the density matrix given in Eq. (|6.10 ), the entropy becomes 



S = 2 I (cosh r]) 2 ln(cosh rj) — (sinh r]) 2 ln(sinh 77) j . (6. 13) 

This expression can be translated into a more familiar form if we use the notation 

tanh?7 = exp ( — — J . (6.14) 



kTJ 

The ratio u/kT is a dimensionless variable. In terms of this variable, the entropy takes the 
form 

s =(^)^RW^T-' n[1 - exp( -" AT)l ' (6 ' 15) 

This is the entropy for a system of harmonic oscillators in thermal equilibrium. Thus, for 
this oscillator system, we can relate our ignorance to the temperature. 
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VII. CONCLUDING REMARKS 



It is interesting to note that Feynman's rest of the universe appears as an increase in 
uncertainty and entropy in the system in which we are interested. In the case of coupled 
oscillators, the entropy allows us to introduce the variable which can be associated with the 
temperature. The density matrix is the pivotal instrument in evaluating the entropy. At 
the same time, the Wigner function is convenient for evaluating the uncertainty product. 
We can see clearly from the Wigner function how the ignorance or the increase in entropy 
increases the uncertainty in measurement. 

The major strength of the coupled oscillator problem is that its classical mechanics is 
known to every physicist. Not too well known is the fact that this simple device has enough 
symmetries to serve as an analog computer for many of the current problems in physics. 
Indeed, this simple system can accommodate the symmetries contained in 0(3,3) which 
is the group of Lorentz transformations applicable to three space-like and three time-like 
dimensions 0. The group 0(3,3) is has many interesting subgroups. Many, if not most, of 
the symmetry groups in physics are subgroups of this 0(3,3) group. 

The authors are grateful to Leehwa Yeh for many helpful comments and criticisms, and 
for pointing out errors in a number of equations in an earlier version of the manuscript. 
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